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Abstract 

In this paper, we propose an inventory model where items are inspected 
through multiple screening processes before delivery to customers. Each screen- 
ing process has independent screening rate and defective percentage. Defective 
items screened out are stored and then returned to supplier. Shortage back- 
ordering are also allowed in the model. Two approaches are used to obtain 
the closed-form optimal order size and the maximum backordering quantity. 
Numerical examples are also provided to demonstrate the use of the model. 

1 Introduction 

When a batch of items arrives from supplier, imperfect quality items do exist in the 
batch. Th erefore, screening process es are required before the items are delivered to 
Salameh fc Jaberl (120001 ) first extended the classical economic order quan- 



customers. 



tity (EOQ) model by integrating a screening process for imperfect quality items into 
th e model. Since then man y researchers proposed exten s ions t o the model developed 
by ISalameh fc Jaberl (120001 ) . iGoyal fc Cardenas-Barron proposed an approxi- 

mation ^o_tJie_jriod^l_by simpjifying the expected reven ue and expected c o st per unit 



time. 



Maddah fc Jaber 



(hoosh pointed out a flaw in ISalameh fc Jaberl (120001) anc 



reformulated the model using the renewal reward theorem. IWahab fc Jaberl (120101 ) 
considered the case when there are different holding costs for the imperfect quality 
and scree ned items. Other e xtensions of this original model can be found in the review 
paper by 



Khan et al. 



(2011) 
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One natural extension to an EOQ model is to consider shortage situation. 



WeeetaL 



(120071 ) develo ped an inventory m odel for items with imperfect quality and shortage 
backordering. Chang &: Hoi ( 2010 ) then revisited their study and provided closed- form 
solutions using the renewal reward theorem. However, in both works they assumed 
that the items sat isfying the backorder are delivered without any screening process. 



Hsu fc Hsul (120121 ) then proposed a model which items are first screened before de- 



livery t o customers . In t his paper, we aim at providing a general version of the 
work of Hs u fc Hsu (bold ) by considering multiple screening processes and shortage 



Pi 



backordering. 

Suppose that after the replenishment of items arrives from the supplier, the items 
have to go through n screening tests before delivering to customers. Denote by Si the 
type i screening process and by Xi the corresponding screening rate of Si. Without 
loss of generality, we suppose 

xi > x 2 > ■ ■ ■ > x n . 

Let pi be the proportion of items in the lot that can not pass Si. Here we assume 
that all Pi are independent of each others in the sense that whether an item can pass 
Si does not depend on the result of other screening processes. Since for any i < j, Si 
will finish before Sj, the proportion of items in y that is screened out after Si is 

Pi, for i = 1; 

£[(l-p fc )pj, for i = 2, . . . ,71. 

k=l 

Then the total proportion of items in y that are classified as defective is given by 

n 
i=l 

The items go through screening processes one by one. Here we assume that the 
defective items are evenly distributed in the lot. Hence, the screening time for the 
Si can be modeled as y/xi. The screening processes are illustrated in Figure 1. In 
Figure 1, the arrows represent the screening rates and the lines represent the items 
undergoing the corresponding screening process. 

The remainder of the paper is organized as follows. In Section 2, we present the 
mathematical model for the problem and derive the expected profit per unit time 
using the renewal reward theorem. In Section 3, we obtain the optimal order size 
and the maximum backordering quantity with an analysis of the model. We then 
give numerical examples to demonstrate the use of our model in Section 4. Finally, 
concluding remarks are given in Section 5 to conclude the paper and address further 
research issues. 
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2 Mathematical model 



In this section, we first give the notations and assumptions used in the development 
of the mod el. Some o f the n otations and assumptions were also used in 



(2012 


) and 


Wee et al. 


(2007) 



Hsu fc Hsu 



Notations: 

y the order size, 

D the demand rate, 

n number of screening processes in a cycle, 

Si type i screening process (i = 1, . . . , n), 

Xi the screening rate of Si, 

c the purchasing cost per unit, 

K the ordering cost per order, 

Pi the defective percentage of Si in y 

Pi the proportion of items in y that screen out after type i screening process 

(i = l,...,n), 

s the selling price of good quality items, 

v the salvage value of defective items (v < s), 

di the unit screening cost of Si, 

B the maximum backordering quantity, 

b the backordering cost per unit per unit time, 

h the holding cost per unit per unit time, 

hd the holding cost for defective items per unit per unit time, 

* the superscript representing optimal value. 
Assumptions: 

1. The demand rate is a known constant. 

2. The lead time of inventory replenishment is assumed to be negligible. 

3. All screening processes and demands proceed simultaneously, but the screening 
rates are greater than the demand rate, > D (i — 1, . . . , n). 

4. The defective items exist in lot size y. 

5. To avoid shortages within screening time, we assume that 

D 

P<1-—. 

•En 

6. The defective items are returned to the supplier when the replenishment items 
arrive. 
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7. Shortage is completely backordered. 



8. A single product is considered. 

The inventory level in a replenishment cycle is illustrated in Figure 2. At the 
beginning of the replenishment cycle, all screening processes proceed simultaneously. 
Therefore, the rate of items to complete all the screening processes depends only on 
the lowest screening process rate x n . After completing all the screening processes, the 
items were first shipped to satisfy the demands. Hence, the inventory level decreases 
at a rate of (1 — p)x n . The remaining items are then used to clear the outstanding 
backorders. Hence, the backordering quantity decreases at a rate of (1— p)x n — D. The 
screening time for Si is y/xi. After the screening process Si, the defective items (piy 
units) are transferred to another inventory warehouse. The inventory level of defective 
items is illustrated in Figure 3. The aim of this paper is to develop a mathematical 
model such that the optimal order size and the maximum backordering quantity can 
be obtained. We first give the length of time intervals in our model. 

[ 1 ~p)y. 



(i) The replenishment cycle length T: T = 

(ii) h = 

(iii) t 2 = T - ti = 

(iv) h 



D 

B 

(1-P)y-B 
D 

B 



(v) t4 — t 2 —t 



(1 - p)x n — D' 

(l-p)y-B B 



D (l-p)x n -D' 

Denote by TR(y) and TC(y,B) the total revenue and the total cost per cycle re- 
spectively. TR(y) is the sum of the total sales of good quality items and the amount 
received from the supplier for the return of the imperfect quality items. One has 

TR{y) = (l-p)ys + pyv. 

TC(y, B) consists of the following costs: 

(i) The purchasing cost: cy; 

(ii) The screening cost: 

n i—1 _ n n—1 n 

^[hx (l-^p^y = [^2^ - Yl dk ) y> 

i=l k=0 i=l i=l k=i+l 
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(iii) The holding cost for good quality items: 
~t\ x D tjx (1 - p)x n 



h x 
= h 
= h 



2 + 2 
(t 2 -t 3 ) 2 D t 2 3 (l-p)x 



+ ^2 PiV x — + 1 3 x t 4 x D 



i=i 

n 



2 2 +J2p*y^ +t ^-^ D . 

t\D tl((l-p)x n -D) ^ y] 



+ 



= h 



2 2 
((l-p)y-B) 2 
2D 



+ 



+ 

i=i 
£? 2 



2((l-pK-L>) 



(iv) The holding cost for defective items: 



i=l i=l 



i=l * 



(v) The shortage cost: 



+ 



2((1 -p)x n -D) 2D 
2£)((l-p)- J D/x n ) 



B 2 ; 



(vi) The ordering cost: K. 

For notational convenience, let 

P 1 = E[(\-p)% P 2 = E\p(\-p)\, P 3 



-, P 



E[p] 



i-E[ P y 4 i-E[ P y 



which depend on p only. Let 



R = 



1 - E[p] 



E 



l-p 



(l-p)-D/x r 



which depends on p, D and x n . Let 



X i 
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which depends on pi, Xi (i = 1, . . . , n) and 



n-l 



A 2 = j2di-J2 E ^{ E * 

i=l i=l k=i+l 



which depends on pi, di (i = 1, . . . , n). Then the expected net profit per cycle is given 
by 

ETP(y,B) 



(1 - E[p])ys + E[p}yv -{cy + A 2 y 
+h 



Pi 2 1 - E[p) n B 2 . 2 
2D y ir lBy+ 2RF-D +My . 



+ h a 



-B 2 + K 



and the expected replenishment cycle length -E[T] 



{} - E[p})y 

D 



By the renewal reward theorem, the expected profit per unit time is 

ETP(y,B) 



ETPU(y,B) 



E[T] 

sD + vDP 4 - (c + A 2 )DP 3 - ( ^ - liB 



+ 



hPi 



+ hA 1 D + h d {P 2 - A 1 D) 



2Ry 

P 3 y + 



KDP 3 

y 



Our aim is to find the optimal order size y and the maximum backordering quality 
B such that the expected profit per unit time ETPU(y, B) is maximized. Since the 
first three terms in (CQ) are independent of B and y, the optimization problem reduces 
to minimizing f(y, B) which is given by 



+ 



2Ry 
hP x 



+ hA x D + h d {P 2 - AiD) 



P 3 y + 



KDP-i 



(2) 



y 



3 Analysis 



In this section, we present two different approaches to minimize f(y,B) in (J2J). 
Approach 1 



We use an approach similar to IChang fc Hoi (120101 ) here. Notice that (j2J) can be 
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rewritten as 



h + b/„ hR \2 
B - - ry 



2Ry 



h + b' 



+ hP 1 + 2hA 1 D + 2h d (P 2 - A X D) 
KDP 3 



h 2 R 



(h + b)P 3 



Pzy 



+ 



y 

Hence, for any fixed y, f(y, B) is minimized when 

hR 

= ~h + l y ' 

Hence the problem reduces to minimizing 

f(y) = [hPi + 2hA 1 D + 2h d {P 2 - A X D) - 
Applying the AM-GM inequality gives 



h 2 R 



(h + b)P 3 



PzV | KDP 3 



f(y) > Pa\ 2KD hP 1 + 2hA 1 D + 2h d {P 2 - A X D) 



h 2 R 



(h + b)P 3 



The equality holds when 



hP x + 2hA 1 D + 2h d (P 2 - A X D) 



which means f(y) is minimized when 



h 2 R 



(h + b)P 3 



P 3 y _ KDP 3 
2 y ' 



2KD 



From (HI), we have 



hP x + 2hA 1 D + 2h d (P 2 - A X D) 



B* = T y 



h 2 R 



(h + b)P 3 



h + b 

Hence the minimum value of f(y, B) is 



f(y*, B*) = P 3 \ 2KD 



hPi + 2hA x D + 2h d (P 2 - A X D) - 



h 2 R 



(h + b)P 3 



and the maximum expected profit per unit time is 
ETPU* = sD + vDP i -(c + A 2 )DP 3 



-P.J2KD 



hPi + 2hA 1 D + 2h d (P 2 - A X D) - 



h 2 R 
{h + b)P 3 



Approach 2 

Differentiating f(y,B) in (j2J) with respect to B and y respectively give 



df_ _ (h + b)B 
dB ~ Ry ' 



df (h + b)B 



2 



+ 



hP x 
C1T 



+ hA\D + h d {P 2 - A 1 D) 



dy 2Ry 2 
The second order partial derivatives are 

d 2 f _h + b 
dB 2 ~ Ry ' 

d 2 f _{h + b)B 2 2KDP 3 
dy 2 Ry 3 y 3 
d 2 f (h + b)B 



KDPi 



2y 2 



dBdy Ry 



2 



Notice that 



d 2 f d 2 f ( d 2 f \ 2 _ 2(h + b)KDP 



dB 2 dy 2 \ dBdy J Ry 4 

d 2 f d 2 f d 2 f ( d 2 f \ 2 

Since we assume that p < 1 — D/x n , we have — — - > and t^t^t^ ~ ( t^tt^ ) > 0- 

oB 2 oB 2 oy z V dBdy / 

This implies that f(y, B) is strictly convex for positive B and y. Hence the unique 

df df 

global minimum for positive B and y can be obtained by solving — — = and — = 0, 

dB dy 

which gives the same results as in (jB]) and (j7|). 

Based on the results above, we then present some special case analyses. 

(i) From ([8]), the defective items are with salvage value vDP^ and the holding cost 
for the defective items is 2P 3 y 'KDh&iPi — AiD). Hence there will be no benefit 
to store the defective items if 



m 

(ii) If the supplier does not accept returns, i.e. v = 0, then there is no storage for 
the defective items {ha = 0). The unit selling price should be higher than 



,2Kh/ . „ hA x \\ 



so that positive profit can be guaranteed, i.e. ETPU* > 0. 

We remark that in the literatures mentioned below, the defective items are not 
stored for return. Hence we consider hd = in the following cases for comparison. 



S 



(iii) When n = 1, the proposed model reduces to the one presented in 
(120121 ). with 



Hsu fc Hsu 



\ 



2KB 



h 



Pi + 2AxB 



hR 



(h + b)P 3 



and 



B* 



hR 

h + b l 



(9) 



(10) 



Hsu fc Hsu 



We re mark that the a bove express i ons ar e simpler than those obtained in \ 
(120121 ) . As stated in iHsu fc Hsul (120121 ). if the defective percentage px follows 
uniform distribution with probability density function 



1 

T 
0. 



for < pi < B; 
otherwise. 



then 







E 



1-pi 



ll-px- D/xtl 



1 D In 

(3xi 



B 2 

l-D/x l 
l-P-D/x] 



In what follows, we give an analysis which is absent in IHsu fc Hsul (120121 ) . The 
equations ([9]) and ( ITU]) can be reduced to 



2KB 



(3 2 BfB h(l - (3/2)' 



(h + b) 



and 



B* 



h{2 - g) 
2(h + b) 



D 



In 



D f 1 - D/xi \ 
. + /3xt 11 \l-p-D/xJ 

1 - D/xt 



-1 



y ■ 



:i2) 



For the case when B is small, by the Taylor series expansion, we can use the 
approximation 



Ill 







.1-B-D/xxJ l-P-D/xx 
By neglecting the terms with (3 2 , we can further reduce and (fT0|) to 



2KD(h + b) 



and 



h[h(l - f3)D/ Xl + b(l-B- DB/xx)} 

h(2-/3)(l-B-D/xx) * 
1 2 (/i + 6)(l-/3) 2/1 ' 
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(13) 



(14) 





Si 


s 2 


s 3 




s 5 


s 6 


S 7 


Pi 


0.01 


0.04 


0.1 


0.01 


0.04 


0.04 


0.1 


Xi (unit/min) 


1 


1 


1 


2 


2 


0.5 


0.5 


di ($/unit) 


0.5 


0.5 


0.5 


1 


0.5 


1 


0.3 



Table 1: Parameters for the screening processes Si (i = 1, . . . , 7). 



Finally, if all products are of good quality, which means the screening rate can be 
set as x\ — > oo, then we have E[(l — p) 2 ] = 1, E[p] = and E 1 



1. Hence y* reduces to the classical EOQ with shortages, 
(iv) When b — > oo, i.e. no shortages are allowed, then 



1 — p — D/xi 



y *-> 



2KD 



+ 2A 1 D) 



which is a generalization of the result obtained in iMaddah Sz Jaberl (120081 ) 

4 Numerical examples 



In this section, we present some numerical examples to demonst rate the use of our 
mo del. We apply the following parameters which are also used in iHsu fc Hsu! (120121 ) 
and lwee et all ( boOT^ . 

Demand rate, D 

Ordering cost, K 

Holding cost, h 

Purchasing cost, c 

Backordering cost, b 

Selling price of good quality items, s 

The salvage value of defective items, v 
Suppose that there are seven types of screening processes, Si (i = 1, . . . , 7). The 
corresponding defective percentage, Pi (i = 1,...,7), are assumed to be uniformly 



50,000 units/year, 

100/cycle, 

$5/unit/year, 

$25/unit, 

$10/unit/year, 

$50/unit, 

$20/unit. 



distributed with respective pdf 




for < pi < A; 
otherwise. 



The parameter the screening rate Xi and the screening cost di for the screening 
process Si (i = 1, ... ,7) are given in Table 1. The optimal y* and B* from and (J7|) 
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y 

B* 

ETPU{y*,B*) 
V* 



Jh 

1624.85 

384.34 

1217432.76 

1630.52 

384.90 



J2 

1638.40 

379.32 

1213159.67 

1662.38 

381.61 



Jh 

1664.90 

368.63 

1204203.81 

1732.15 

374.58 



Jh 

1679.81 

477.24 

1192509.48 

1682.95 

477.72 



_5s 

1699.16 

474.22 

1213382.37 

1712.65 

476.31 



So 



St 



1534.16 1542.35 

209.17 194.28 

1187226.30 1214227.78 

1573.85 1651.62 

208.45 191.32 



ETPU{yl,Bl) 1217432.72 1213158.94 1204198.33 1192509.47 1213382.17 1187223.69 1214208.42 



Table 2: Numerical results for the screening processes Si 



7). 



and the corresponding expected profit per unit time for each of the screening process 
are given in Table 2. We also give the approximated y* and B* from ( fl3|) and ( fl4|) 
to demonstrate the effectiveness of the approximations. The results show that the 
approximations are good as the errors of the expected profit per unit time are all less 
than 0.01%. The results of S\, S2 and 63 show that when (3 increases, y* increases 
but B* and ETPU* decreases. The results of S2 and S5 show that when x increases, 
all y*, B* and ETPU* increases. The results of S3 and S7 show that when both /3 
and d increase, both y* and B* increase but ETPU* decreases. 

We then provide the results for different combinations of the screening processes. 
First we consider the case when there are two screening processes in the cycle. The 
results are given in Table 3. We then compare these results with those in Table 2. 
We remark that for the case of two screening processes, the optimal order size y* is 
between those in the cases of one screening process. However, when there are two 
screening processes in the cycle the maximum backordering quantity B* is less than 
both B* in the cases of one screening process. 

Next we consider the case when there are three screening processes in the cycle. 
The results are given in Table 4. We observe that, for example in the case of Si + 
S4 + S e , the optimal order size y* is between those in the cases of Si + £4 and S4 + Sg; 
and it is also close to that of 5*6, which has smallest Xi among Si, S4 and Sq. The 
maximum backordering quantity B* is less than those in the cases of Si + 5*4 and 
S4 + Sq] and it is also close to that of Sq, which has smallest Xi among Si, S4 and Sq. 
We remark that similar patterns can also be observed in other cases. 

5 Concluding Remarks 

In this paper, an inventory model for items with imperfect quality are developed. 
The items are inspected through multiple screening processes and shortage are back- 
ordered. The optimal order size and maximum backordering quantity are obtained 
by two approaches. We then provide an analysis of the model. Numerical examples 
are also provided to demonstrate the use of the proposed model. The current study 
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Si +S 4 


S-2 + S4 


53 + 54 


S1 + S5 


S 2 + S B 


S3 + S5 


y* 


1630.93 


1644.44 


1670.86 


1649.22 


1662.63 


1688.81 


B* 


383.06 


378.01 


367.26 


379.11 


373.97 


363.01 


ETPU{y*,B") 


1165658.46 


1160592.66 


1149976.27 


1186633.71 


1181888.2 


1171942.72 






Si+S 6 


s 2 + s 6 


s 3 + s 6 


Si + S 7 


S2 + S7 


S3 + S7 


y* 


1538.28 


1550.54 


1574.43 


1546.29 


1557.99 


1580.69 


B* 


207.19 


201.04 


187.85 


192.23 


185.89 


172.20 


ETPU{y*,B") 


1160290.47 


1155925.26 


1146776.05 


1186440.53 


1181934.97 


1172491.58 



Table 3: Numerical results for different combinations of the screening processes (n = 
2). 





Si + S4 + S(, 


S2 + S4 + S$ 


S3 + S4 + Se 


Si + S 5 + S 6 


s 2 + s 5 + 


s 6 


S3 + S5 + Se 


y* 


1543.69 


1555.9 


1579.7 


1559.95 


1572.02 




1595.49 


B* 


205.36 


199.17 


185.89 


199.69 


193.37 




179.78 


ETPU{y*,B") 


1107457.62 


1102283.46 


1091439.66 


1128854.96 


1124007.4 




1113847.66 






Si + s 4 + s 7 


S2 + S4 + s 7 


S3 + S 4 + S 7 


Si + s 5 + s 7 


S 2 + S s + 


s 7 


S3 + S5 + S 7 


y* 


1551.53 


1563.17 


1585.73 


1567.27 


1578.71 




1600.83 


B* 


190.33 


183.94 


170.15 


184.43 


177.89 




163.74 


ETPU(y*,B*) 


1131938.47 


1126598.34 


1115406.7 


1154009.02 


1149005.6 




1138519.14 



Table 4: Numerical results for different combinations of the screening processes (n = 



can be extended in several ways. In this study, all screening processes are of 100% 
accurancy. The screening rates are also assumed to be greater than the demand rate. 
These assumptions may be relaxed in future research. Another direction is to consider 
stochastic demand rates. 
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A X2 A± x 3 



S\ S2 S3 S n 
Figure 1: The screening processes. 



inventory level 




time 



T 



Figure 2: The inventory level in a replenishment cycle. 
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inventory level 



PnV 

pry | i ! 

y_ y_ y_ t 

Xi x 2 x n 
Figure 3: The inventory level of defective items. 



15 



